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Abstract  -  This  paper  presents  a  coupled,  joint 
probabilistic  data  association  (JPDA)  algorithm  for 
multi-target  tracking  using  a  modified  version  of  the 
standard  measurement-to-track  assignment  model. 
The  mutually  exclusive  nature  of  standard  JPDA 
association  events  precludes  any  measurement  be¬ 
ing  associated  with  more  than  one  target  in  a  given 
event.  This  constraint  is  relaxed  here  to  allow  a  mea¬ 
surement  to  be  assigned  to  multiple  targets.  All  other 
JPDA  assumptions  are  retained  (i.e.,  no  measure¬ 
ment  can  be  simultaneously  associated  with  target 
and  clutter,  and  each  track  can  claim  at  most  one 
measurement).  The  computational  requirements  of 
the  resulting  algorithm  grow  linearly  with  the  num¬ 
ber  of  tracks.  The  recursive  estimators  for  the  cou¬ 
pled  track  means  and  covariance  are  derived  and  pre¬ 
sented. 
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1  Introduction 

The  joint  probabilistic  data  association  (JPDA)  al¬ 
gorithm  [1-4]  for  multi-target  tracking  in  clutter  has 
many  desirable  properties.  Prominent  among  these 
properties  is  the  fact  that  the  track  estimates  are  ob¬ 
tained  from  the  posterior  distribution  of  the  states, 
thereby  providing  a  more  realistic  assessment  of  the 
uncertainty  in  the  multi-target  tracking  problem.  The 
single  target  version,  the  probabilistic  data  associa¬ 
tion  filter  (PDAF)  [5],  is  a  consistent  estimator  [2,  3,  6] 
that  works  well  in  many  different  applications.  In  fact, 
PDAF  is  regarded  by  many  to  be  the  baseline  for  track¬ 
ing  a  target  in  clutter. 

The  primary  drawback  of  JPDA  is  the  explosive 
growth  in  association  events  that  must  be  considered 
as  the  number  of  targets  and  detected  measurements 
increase  to  even  moderate  levels  («10).  Approxima¬ 
tions  to  JPDA  that  avoid  this  explosive  growth  are 
discussed  in  references  [7-11].  This  paper  presents  a 


modified  JPDA  algorithm  whose  number  of  association 
events  grow  linearly  with  the  number  of  targets,  as  op¬ 
posed  to  the  exponential  growth  in  JPDA.  Recently,  a 
new  linear  complexity  multi-target  tracking  algorithm 
based  on  the  integrated  probabilistic  data  association 
algorithm  (IPDA)  [12]  and  joint  IPDA  (JIPDA)  [13] 
was  described  in  references  [14-16].  The  algorithm  in 
[14-16]  is  fundamentally  different  from  the  algorithm 
presented  here  because  the  assignment  assumptions  are 
different.  Performance  comparisons  between  the  two 
algorithms  are  currently  being  initiated  and  will  be 
presented  at  a  later  time. 

The  modified  JPDA  algorithm  is  developed  in  the 
next  several  sections.  The  assignment  model  is  pre¬ 
sented  in  Sec.  2.  The  filtering  algorithm  that  results 
from  this  assignment  model  is  discussed  at  a  high  level 
of  abstraction  in  Sec.  3,  and  these  results  are  special¬ 
ized  to  linear  Gaussian  targets  and  uniform  clutter  in 
Sec.  4.  An  example  containing  two  crossing  targets  is 
then  examined  in  Sec.  5,  followed  by  a  short  summary 
and  conclusions  in  Sec.  6. 

2  Assignment  Model 

In  the  following  derivation  and  discussion,  targets  are 
assumed  to  be  close  enough  together  so  that  a  multi¬ 
target  tracking  procedure  is  required  (e.g.  the  target 
gates  overlap).  For  widely  separated  targets  or  a  sin¬ 
gle  target  in  clutter,  the  PDAF  is  employed.  For  the 
closely  spaced  targets,  the  JPDA  algorithm  is  applied, 
but  with  a  relaxed  set  of  assignment  assumptions.  In 
particular,  whereas  JPDA  does  not  permit  a  measure¬ 
ment  to  be  associated  with  multiple  targets  in  a  sin¬ 
gle  event,  the  new  model  allows  measurements  to  be 
assigned  to  multiple  targets  (but  not  target  and  clut¬ 
ter)  by  assuming  that  the  measurement-to-track  as¬ 
signments  are  statistically  independent  between  tar¬ 
gets.  A  mutually  exclusive  and  exhaustive  set  of  as¬ 
sociation  events  is  considered  for  each  target. 

To  establish  notation,  the  event  in  which  the  /cth 
measurement  is  assigned  to  the  jth  target  track  dur- 
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ing  the  nth  scan  interval,  and  in  which  all  other  mea¬ 
surements  are  assigned  to  clutter,  is  denoted  by  Ojkn- 
The  vector  represents  the  collection  of  association 
events  for  the  jth  target  track.  Consequently, 

“  LJ 

k=0 

where  p|  6 jin  =  0  for  all  fc  yf  ?, 

Prob  ^9ikn  Pi  =  Prob  {e,kn}  Prob  {6»y/„}  (2) 

for  t  yf  j  and  Kn  represents  the  number  of  measure¬ 
ments  in  scan  n.  The  association  events  for  each  tar¬ 
get  make  no  explicit  allowance  for  the  existence  of 
other  targets  (unlike  the  methods  described  in  [14-16]). 
Multi-target  behavior  arises  due  to  the  mixing  that  re¬ 
sults  from  the  independence  of  the  between-target  as¬ 
sociation  events. 

The  measurements 

Z„  =  (3) 

{Kn  >  0)  in  scan  n  are  statistically  independent  when 
conditioned  on  the  target  states,  Xj„  (for  j  =  1, . . . ,  J), 
and  the  association  events,  9jn-  Furthermore,  the  tar¬ 
get  states  are  independent  of  the  association  events. 
Stacking  the  state  vectors  from  all  targets  into  a  single 
vector  defines  the  combined  state  vector 

X„  =  [xL...,xL]^  (4) 

whose  evolution  from  scan  to  scan  is  governed  by  a  set 
of  independent  first  order  Markov  processes,  such  that 

,7 

p(X„|X„_i)  =  p(xy„|xj„_i).  (5) 

j=i 

Here,  J  denotes  the  number  of  targets  in  the  obser¬ 
vation  space  with  overlapping  target  gates.  While 
the  measurements  are  conditionally  independent  when 
given  the  combined  state  vector,  the  individual  target 
state  estimates  are  not  necessarily  independent  when 
target  gates  overlap.  Dependent  target  state  estimates 
are  said  to  be  coupled,  and  they  are  estimated  jointly 

[3]. 

The  posterior  event  probabilities  for  the  modified  al¬ 
gorithm  are  obtained  via  the  usual  application  of  Bayes 
theorem.  The  measurement  likelihood  for  event  9jkn, 
in  which  the  fcth  measurement  goes  to  the  jth  target, 
is  given  by 

P{'^n\^jm  j)  ~  P{'^kn\^jn)  Pc{^rn) 

r—1 
r^k 


where  Pc{-)  denotes  the  likelihood  of  a  clutter  mea¬ 
surement.  For  the  event  that  assigns  all  measurements 
to  clutter  (i.e.,  fc  =  0),  the  measurement  likelihood  is 
given  by 

K„ 

P{'^n\9j0mj)  —  Pci'^kn)-  (7) 

fc=l 

The  joint  likelihood  of  the  measurements  conditioned 
on  the  state  for  target  j  is  obtained  by  marginalizing 
over  the  association  events  for  target  j  to  obtain 

p(Zn|Xj>j,J)  =  ^  ]  p{’Zin\^jn:  9jkn'!  j)  P^ob  ■ 

k=0 

(8) 

This  expression  is  the  conditional  measurement  likeli¬ 
hood  for  the  (single-target)  PDAF.  Similarly,  the  like¬ 
lihood  of  the  measurements  in  scan  n,  conditioned  on 
all  the  target  states,  is  obtained  by  marginalizing  over 
the  targets  as 

J 

p(Z„|X„)  =  ^  Probjj}  p(Z„|xj„,j)  (9) 

where  Prob{j}  represents  the  prior  probability  of  tar¬ 
get  j. 

3  Posterior  Estimates 

The  estimates  of  the  target  states  that  minimize  the 
posterior  mean  square  error  and  the  corresponding  co- 
variance  matrices  at  each  scan  are  obtained  by  tak¬ 
ing  expectations  with  respect  to  the  conditional  target 
state  densities.  The  process  of  computing  these  expec¬ 
tations  is  broken  into  two  steps  consisting  of  a  time 
update  and  a  measurement  update  as  described  in  [17] 
and  [18].  Just  enough  of  the  theory  is  repeated  in  this 
section  to  introduce  needed  notation. 

The  starting  point  for  developing  state  estimators  is 
the  joint  likelihood  function 

P(Zn5  1 

=  p(Z„|X„)p(X„,X„_i|Z„_i) 

=  p(Z„|X„)p(X„|X„_i)p(X„_i|Z„_i)  (10) 

where  Z„_i  represents  the  collection  of  all  the  mea¬ 
surements  from  the  starting  scan  through  and  includ¬ 
ing  scan  n  —  1.  Given  the  (prior)  PDF  of  the  previous 
state  estimates,  p(X„_i|Z„_i),  the  time  update  forms 

p(X„|Z„_i)  =  J  dXn-lp{Xn,Xn-l\Zn-l).  (11) 


(6) 


The  measurement  update  applies  Bayes  rule  to  obtain 

p(X„|Z„)  =  p(X„|Z„,Z„_i) 

^  p(Z„|X„)p(X„|Z„_i) 

p(Z„|Z„_i)  ■  ^  ^ 

The  denominator  term  p(Z„|Z„_i)  is  the  combined  in¬ 
novation  likelihood: 

p(Z„|Z„_i)=  y'dX„p(Z„,X„|Z„_i).  (13) 

In  practical  applications,  p(Z„|Z„_i)  is  based  on  a 
truncated  density  function  such  that  only  measure¬ 
ments  near  the  predicted  target  measurement  are  con¬ 
sidered  by  the  algorithm.  This  procedure  is  referred  to 
as  gating  or  measurement  validation  [2-4, 19]. 

Given  the  posterior  state  density,  the  state  is  esti¬ 
mated  as 


X 


7i\n 


dX„  X„p(X„|Z„). 


(14) 


distribution  of  is  obtained  by  convolving  the  bi¬ 
nomial  distribution  for  the  number  of  target  measure¬ 
ments  with  the  Poisson  distribution  for  the  number  of 
clutter  measurements,  which  gives 


Prob{X„}  = 

e--f 

1=0 


(17) 


where  J'  =  min(itr„,  J).  Using  Eq.  (17),  and  following 
the  derivation  in  Appendix  D.4  in  [2],  the  prior  prob¬ 
ability  for  the  event  in  which  the  ki\i  measurement  is 
assigned  to  the  jth  target  (and  all  other  measurements 
to  clutter)  is  given  by 


Proh{6jkn}  = 

L-i  f/i(^n)/^n  forl<fc<K„ 

\fo{Kn)  ioT  k  =  0 


Prob{Ar„}' 


(18) 


where 


Defining  the  outer-product  function  ^(X)  =  XX"*", 
the  posterior  covariance  matrix  of  the  combined  tar¬ 
get  state  vector  is  given  by 

7^„l„  =  J  dX„  5(X„  -  X„|„)  p(X„|Z„).  (15) 

4  Linear  Gaussian  Targets  in 
Poisson  Clutter 

This  section  develops  the  time  and  measurement  up¬ 
dates  for  linear-Gaussian  targets  and  Poisson  clutter. 
Target  measurements  are  detected  with  a  probability 
of  Pd,  and  clutter  measurements  are  uniformly  distrib¬ 
uted  in  the  observation  space.  The  modified  JPDA 
algorithm  considers  an  observation  space  that  is  the 
union  of  the  overlapping  target  gates.  Denoting  by  V 
the  volume  in  the  union  of  the  gates,  the  clutter  density 
in  this  observation  region  is 

Pc(zfc„)=V-i.  (16) 

The  number  of  clutter  measurements  is  modeled  by  a 
Poisson  distribution  with  parameter  AV,  where  A  is 
the  average  number  per  unit  volume  (spatial  density). 
The  number  of  measurements,  Ar„,  in  a  scan  is  equal  to 
the  sum  of  the  number  of  target  measurements  and  the 
number  of  clutter  measurements  which  are  statistically 
independent. 

The  number  of  target  measurements  in  any  given 
scan  is  between  zero  (no  target  produces  a  measure¬ 
ment)  and  J  (all  targets  produce  a  measurement)  and 
has  a  binomial  distribution  with  parameter  Pd-  The 


fi{Ku)  = 


,-Ay' 


1=1 


^(f)p'(l-P,)-^' 


Kn-l 


(7f„-/)! 


(19) 


and 

/o(i^„)=e-^^(l-Prf)^'^^p^.  (20) 

When  only  one  target  is  present  Eq.  (18)  reduces  to  the 
event  prior  probability  given  for  the  PDAF  in  [2,  3] . 

The  posterior  event  probabilities  are  obtained  by 
including  the  appropriate  innovation  likelihood  terms. 
Under  linear  Gaussian  assumptions,  the  likelihood  of 
a  target  measurement  equals 


P(2fcn|^jn)  —  A/" 

~  A7  (21) 

where  Af  (x;  fi,  a)  denotes  the  Gaussian  density  with 
mean  /i  and  covariance  a  and 

B,„  =  [0  •  •  •  0  B,„  0  •  •  •  0] .  (22) 


The  density  function  for  a  target’s  current  state,  given 
the  previous  state,  is  defined  as 

l)  —  Af  {^jn,  .  (23) 

The  joint  density  of  all  target  states  at  scan  n,  condi¬ 
tioned  on  all  the  data  up  to  and  including  scan  n,  is 
also  Gaussian  and  is  given  by 

p(X„|Z„)=A7(x„;X„|„,7^„I„).  (24) 

Because  the  covariance  matrix  7^„|„  is  not  block  diag¬ 
onal  (as  will  be  shown  subsequently)  the  target  state 
estimates  are  coupled. 


4.1  Time  Update 

For  linear  Gaussian  targets,  Eq.  (5)  is  equivalent  to 

p(X„|X„_i)=AA(X„;.4„X„_i,Q„)  (25) 

where  An  =  diag  {Aj„}  is  the  block  diagonal  ma¬ 
trix  of  the  target  state-feedback  matrices  and  Q„  = 
diag  {Qjn}-  The  time  update  therefore  amounts 
to  marginalizing  the  previous  state  from  the  joint 
Gaussian  density 

p(X„,X„_i|Z„_i) 

=  p(X„|X„_i)p(X„_i|Z„_i) 

=  A/'(X„;  A„X„_i,  Q„) 

X  Af  ^X.n-i;Xn-ij„-i,Vn-lln-l'j  ■  (26) 

For  this  type  of  Gaussian  product,  the  dependency  re¬ 
versal  from  Bayes  theorem  is  easily  achieved  by  apply¬ 
ing  the  Gaussian  refactorization  lemma  (GRL)  [20], 
which  directly  yields  the  expression 

p(X„,X„_i|Z„_i) 

=  p(X„|Z„_i)p(X„_i|X„,Z„_i) 

~  Af  ^X^;  Af  (X^_i^  TJ7^,  A^) 


(27) 

where 

^n\n—l  '^n^n—l\n—l  (28) 

l^n|n-l  =  Qn  +  AnVn-l\n-l  An  (29) 

Un  =  W„  X„_]^|„_]^  -|-  H„  X„  (30) 

A„  =  'WnVn-l\n-l  (31) 

W„  =  I-H„A„  (32) 

and 

^n='Pn-l\n-l^n'Pn\n-l-  (^^) 


With  this  dependency  reversal,  integrating  the  right 
side  of  Eq.  (27)  with  respect  to  X„_i  is  trivial,  giving 

p(X„|Z„_i)  =  Af  ^X„;  X„|„_i,  (34) 

which  is  the  familiar  predicted  state  density.  Since  the 
covariance  matrix,  of  the  combined  state  es¬ 

timate  X„|„_i  is  not,  in  general,  block  diagonal,  the 
predicted  target  state  densities  are  not  statistically  in¬ 
dependent. 


4.2  Innovation  Likelihood 

The  fact  that  the  predicted  target  states  are  coupled 
complicates  the  derivation  of  p(Z„|Z„_i)  to  some  de¬ 
gree.  To  address  this  problem,  the  combined  state  out¬ 
put  matrix  was  introduced  in  Eq.  (21).  By  using 
Bj„X„  instead  of  Bj„Xj„,  the  GRL  can  be  applied  di¬ 
rectly  to  (34)  such  that  marginalization  with  respect 
to  the  combined  target  state  is  straightforward. 

Since  the  clutter  terms  in  Eq.  (6)  do  not  de¬ 
pend  on  the  state  vector  for  target  j,  the  product 
p(Z„|xj„,6»jfc«)p(X„|Z„_i)  is  marginalized  by  apply¬ 
ing  the  GRL  to  the  two  Gaussian  components,  inte¬ 
grating  X„  from  the  result,  and  then  multiplying  by 
the  clutter  terms  after  the  fact.  The  Gaussian  compo¬ 
nents  are  re-factorized  as 


Af  Af  ^^^n^^^n\n—lif^n\n—l 


= 

•N"  ZjTT,,  Sjyj)  A/*  {yLn]~^jk,n\n')'^ jk,n\n) 

(35) 

where 

^jn  ^jn^n|n— 1 

(36) 

^  X 

Sjn  —  R-jn  Bjn^n|n— 1  "^jn 

(37) 

‘^jk,n\n  ^n|n— 1  ^jn  '^kn 

(38) 

jk,n\n  ^jn'^n\n—l 

(39) 

^jn  I  ^jn  ^jn 

(40) 

and 

Grjn  ^n|n— 1  ^jn  ' 

(41) 

Integrating  Eq.  (35)  with  respect  to  X„  gives 

P{'^kn\^jn^  —  P{'^kn\^jn)  p(Xyj  |X^|,^_2) 

=  Af  {^kn'i  ^jn^  ^jn)  •  (d2) 

Inclusion  of  the  clutter  terms  then  gives  conditional 
likelihood  as 

PC^nlf^jknj  ^n—l)  =  V  ”  Af  (Zfc„;  Zj„,  Sj„)  (43) 

for  1  <  fc  <  Kn-  For  k  =  0,  where  all  measurements 
are  presumed  to  be  clutter,  this  expression  reduces  to 

p{Zn\9okn,Zn-l)  =  (44) 

The  overall  innovation  likelihood  is  obtained  by  nested 
application  of  the  total  probability  theorem  as 

,7 

p(Z„|Z„_i)  =  ^Prob{j} 
i=i 

Kn 

X  ^  Piob  {Ojkn}  p{Zn\0jkn,Zn-l).  (45) 


4.3  Measurement  Update 

Substituting  equations  (34)  and  (45)  into  Eq.  (12) 
yields  the  posterior  density  of  the  combined  target 
state  as 


p(X„|Z„_i) 


P(^n  |X^)  Af  ,  7^n|n— 

p(Z„|Z„_i) 


(46) 

The  state  estimates  are  obtained  by  taking  expecta¬ 
tions  with  respect  to  this  conditional  density.  In  com¬ 
puting  the  conditional  mean,  there  are  two  types  of 
expectations  that  need  to  be  evaluated.  For  each 
measurement-to-target  association  hypothesis,  the  ex¬ 
pected  value  of  X„  is  given  by 


^jk,n\n  ^(^711-^71—1) 

X  f  cDi-n  Xji  p(Z^|Xjy,,  p{^X.n\^n\n—l^  •  (^^f) 


where  XjQ  „|„  =  X„|„_i.  The  posterior  mean  of  the 
combined  target  state  is  therefore 

X,7|77  =  EE 

j  =  l  k=0 

J  K„ 

~  ^7i|77— 1  T  ^  ^  ^  ^  '^jkn  i^nk  ^jn)  (48) 

i=l  k=l 


where 


Finally,  substituting  Eqs.  (39)-(41)  for  7^„|„  and  sim¬ 
plifying  gives 


—  (I  ~  G„)7^„|„_1 

J  iC„ 

T  ^  ^  ^  '4^jkn  9(^^jk,n\n  ^7i|77^  (^^) 

j=l  k=0 

where 

J 

i=i 

An  alternative  method  of  computing  X„|„  and  Vn\n 
is  obtained  by  breaking  equations  (48)  and  (52)  into 
two  convex  sums  instead  of  one.  The  first  convex  sum 
is  over  measurements  and  is  essentially  a  PDAF  step. 
The  second  convex  sum  mixes  the  PDAF  estimates. 
Approximating  the  weights  used  in  this  mixing  be¬ 
tween  tracks  (second  convex  sum)  by  the  prior  target 
probabilities  yields  a  far  more  stable  algorithm.  To 
begin  the  PDAF  step,  let 


'  K„ 


E/ 


^  jn  ^jri' 


(53) 


Kk=l 


Pjn  =  m  ^n—1  )  Prob{6»jfc77}  (54) 

fc=0 


and 


Pjkn  —  l^jn  -2^71— l)  Fl’Ob  .  (55) 

Then,  the  contribution  to  the  combined  state  vector 
from  target  j  equals 


Probjj}  Prob{6»jfe„}  p(Z„|6>jfc„,Z„_i) 

The  combined  state  covariance  is  obtained  following 
similar  developments.  Recalling  the  definition  ^(X)  = 
XX"*"  and  noting  that 

^^X^^  Xj^  +Xj7. 

the  covariance  matrix  for  each  measurement-to-target 
association  event  is  defined  by 

'^jk,n\n  T  9  (^^^jk,n\n  X^|^^ 

—  P{'^n\-^n—l)  j  dX^^  ^  ^X^ 

^  P(^n|^jn7  ^jkn)  ^(Xt^  |Xj^|yj_  ]^ )  .  (50) 
where  Pjo,n|n  =  d^n|n-i-  Combining  terms  then  yields 


Krt 


'^j,n\n  —  ^  ^  i^jkn^jk. 


n\n 


A;=0 


Kr, 


—  Xyj|^_]^  -h  Gjn  E  Pjkn  {^kn  ^n\n—l^  (56) 


k=l 


which  is  the  same  as  the  PDAF’s  target  state  estimate 
except  that,  in  this  case,  it  is  the  contribution  of  target 
j  to  the  combined  state  vector.  To  compute  Vn\n,  the 
average  distances  from  Xj^  „|„  to  Xj  „|„  and  Xj  „|„  to 
are  required: 

Kn 

~^j,n\n  ^  ^  Pjkn  (^^jk,n\n  Xj  (57) 

fc=l 

and 

X'jj7i|77— 1  =  PjOn  ^Xj  X,^|,.j_x)  .  (58) 

The  covariance  contribution  from  target  j  is  given  by 


'Pn\n  = 

J  Kr, 

E  E 

j  —  1  k—0 


'^j,n\n  —  PjOn  'Pn\n—1  T  9  Xj 

7^jfc,7i|77  +  5^Xjfc  —  X„|„^  .  (51)  +  f^jkn  'Pjk,n\n  +  9  (^jk,n\n  ~^j,n\n^  (59) 


which  is  the  same  form  as  the  PDAF  covariance  esti¬ 
mate. 

To  form  the  single  combined  state  estimate  and  co- 
variance,  we  introduce  the  posterior  probability  of  tar¬ 
get  j,  given  by 


Probjj} 

p{Zn\Zn-l) 


and  observe  that 


j 

=  1- 
i=i 


(60) 


(61) 


Because  ujjn  is  the  posterior  probability  of  target  j, 
and  because  p(Z„|Z„_i)  is  the  average  value  of  the 
/3j„’s  (see  equations  (45)  and  (54)),  it  follows  that 


Update  Rate 

1  Sec 

Number  of  Updates 

65 

Target  1  Course 

88  Deg 

Target  2  Course 

92  Deg 

Target  1  Init  Pos 

(0,  1200)  m 

Target  2  Init  Pos 

(0,  2000  )m 

Target  1  Speed 

500  m/s 

Target  2  Speed 

500  m/s 

Process  Noise  Sigma 

0.01 

Measurement  Noise  Sigma 

75 

A 

1  per  km^ 

Pd 

0.99 

Pg 

0.99 

Table  1:  Crossing  tracks  scenario  parameters 


Ujn  ~  Prob  {j}  . 

Using  this  approximation  for  yields  more  stable 
estimates  because  it  prevents  one  or  more  targets  from 
dominating  the  mixing  across  targets,  which  can  cause 
the  covariance  estimates  for  the  remaining  targets  to 
blow  up  over  time. 

Using  Eq.  (56),  the  estimate  of  the  combined  target 
state  vector  from  Eq.  (48)  is  given  by 

J 

^n|n  ^  (^^) 

i=i 

Similarly,  Eq.  (52)  now  may  be  expressed  as 


MJPDA 


(62) 


J  r 

n\n  ^  ^  ^jn  '^j,n\n  T  9  ^n|n^ 


+  Y,- 


j,n\n  ^n|T 


j,n\n  -^n\n  j  ^  j,n\r, 
T 


(64) 


5  Example 

A  two  target  crossing  scenario  was  used  to  compare  the 
performance  of  the  Modified  JPDA  algorithm  against 
single  target  PDAF  and  a  multiple  target  JPDA.  Per¬ 
formance  criteria  include  the  number  of  lost  tracks,  the 
Normalized  Estimation  Error  Squared  (NEES),  posi¬ 
tional  error,  and  algorithm  complexity.  A  lost  track  is 
defined  as  any  track  whose  NEES  value  exceeds  twenty 
during  any  update.  One  thousand  Monte  Carlo  trials 
were  used  to  generate  the  results.  The  scenario  pa¬ 
rameters  used  for  the  simulation  are  given  in  table  1, 


Figure  1:  Sample  scenario  run  with  MJPDA  results 


Figure  2:  RMS  positional  error  curves  for  PDAF, 
JPDA  and  MJPDA 


Alg. 

%  Lost 
Tracks 

NEES 

Rel.  Pos. 
Error 

Complexity 

PDAF 

18.8 

Consistent 

1 

Linear 

JPDA 

8.0 

Consistent 

1 

Exponential 

MJPDA 

13.0 

Mostly  consistent 

3 

Linear 

Table  2:  Performance  summary 


Figure  3:  NEES  curves  for  PDAF,  JPDA  and  MJPDA 


and  an  example  showing  results  for  MJPDA  appears 
in  figure  1.  The  results  from  the  runs  are  summarized 
in  table  2  and  figures  2  and  3. 

This  example  indicates  that  the  modified  JPDA  al¬ 
gorithm’s  lost  track  performance  is  better  than  that  of 
independent  PDAFs,  but  is  degraded  relative  to  JPDA. 
Furthermore,  its  covariance  estimates  are  too  large 
during  the  cross,  making  the  NEES  measure  too  small 
in  that  region,  and  the  position  errors  are  larger  than 
PDAF  and  JPDA.  The  degraded  performance  charac¬ 
teristics  are  believed  to  result  from  splitting  all  the 
measurements  in  the  union  of  two  intersecting  gates 
between  the  two  targets.  As  a  result,  the  measure¬ 
ments  have  a  reduced  impact  on  the  updates  to  the 
target  states  in  Eqs.  (48)  and  (63)  and  the  target  state 
covariance  matrices  in  Eqs.  (52)  and  (64).  The  pre¬ 
dicted  state  and  covariance  are  therefore  more  influ¬ 
ential,  causing  the  state  covariance  to  grow  too  large, 
which,  in  turn,  causes  the  NEES  plot  for  MJPDA  to 
fall  below  the  lower  threshold  in  figure  3.  This  prob¬ 
lem  may  be  addressed  by  splitting  between  targets  only 
those  measurements  that  fall  within  the  intersection 
(as  opposed  to  the  union)  of  the  target  gates.  An  in¬ 
vestigation  of  that  approach  is  ongoing. 


6  Summary 

This  paper  presents  a  modified,  coupled  JPDA  algo¬ 
rithm  that  achieves  linear  computational  complexity 
by  relaxing  the  requirement  that  a  measurement  can 
only  be  assigned  to  one  target.  Instead,  this  new  algo¬ 
rithm  assumes  that  the  measurement  to  target  assign¬ 
ment  events  for  target  are  statistically  independent,  al¬ 
lowing  a  given  measurement  to  be  assigned  to  multiple 
targets.  In  this  single  target  case,  the  new  algorithm 
reduces  to  the  PDAF. 

The  intent  of  this  work  is  to  develop  a  consistent 
estimator  whose  performance  is  near  to  that  of  JPDA 
but  whose  computational  growth  is  linear  in  the  num¬ 
ber  of  targets  and  measurements  (instead  of  exponen¬ 
tial).  Future  work  will  include  an  examination  of  the 
MJPDA  algorithm  in  which  only  those  measurements 
in  the  intersection  of  gates  is  split  between  targets  (as 
discussed  in  the  previous  section)  and  a  performance 
comparison  between  MJPDA  and  the  IPDA-based  al¬ 
gorithm  [14-16]. 
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